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Abstract 

Relativistic systems of particles interacting pairwise at a distance (interactions not 
mediated by fields) in flat spacetime are studied. It is assumed that the interactions 
propagate at the speed of light in vacuum and that all masses are scalars under Poincare 
transformations. The action functional of the theory depends on multiple times (the 
proper times of the particles). In the static limit, the theory has three components: a 
linearly rising potential, a Coulomb-like interaction and a dynamical component to the 
Poincare invariant mass. In this Letter we obtain explicitly, to terms of second order, 
the Lagrangian and the Hamiltonian with all the dynamical variables depending on a 
single time. Approximate solutions of the relativistic two-body problem are presented. 
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After the discovery of the action-at-a-distance formulation of electrodynamics [T] - [7J, 
several relativistic non-instantaneous action-at-a-distance theories have been investigated [8] 
- |21j . Instantaneous action-at-a-distance formulations have been studied using a variety of 
approaches (22] - [30] . 

In a recent paper [31] the relativistic action-at-a-distance approach of Wheeler and Feyn- 
man [1] - [7J was extended in order to explore what other types of interparticle interactions 
are allowed in special relativity. It was assumed that the interactions propagate at the speed 
of light and that the theory can be described by an action that does not depend on the four- 
vector accelerations or on higher derivatives. A general action functional depending only on 
the four-vector velocities and relative positions of the particles was investigated. Assuming 
that the Poincare invariant parameters that label successive events along the world lines of 
the particles can be identified with the particles' proper times and that all masses are scalars, 
the most general form of the interaction terms in the action was determined explicitly [31] . 

The theory presented in [31] is described by an action which depends on multiple times. 
Quantum mechanics, on the other hand, is based on the hamiltonian formalism (or the Dirac 
hamiltonian formalism for constrained systems [32] - [35]), which is built on the assumption 
that one can choose a single time variable to describe the evolution of the whole system. 
The lagrangian (or covariant) quantization procedure [36] - [38] also assumes a single time 
variable. 

One may attempt to transform a theory with multiple time variables into an equivalent 
theory with just one time variable, in a given inertial reference frame, using Taylor series 
expansions involving the particles' present motions |39j. The problem of finding a single 
time variable to describe a relativistic system of interacting particles has been investigated 
by several authors [10] - [12]. 

In this Letter we present a reformulation of the theory given in [31] as an action-at- 
distance theory with a single time variable. We use series expansions up to terms of second 
order (^). 

We obtain the approximate Lagrangian for this theory. Approximate expressions for the 
total energy, total linear momentum and total angular momentum are presented. We also 
obtain the Hamiltonian to terms of second order. 

We present a formula for the Hamiltonian of two interacting particles in the center of 
momentum reference frame. Using this result, direct predictions for the energy spectrum 
can be made. 

We also present approximate circular orbits solutions of the relativistic two-body problem. 
It is well known that an isolated system of N particles interacting electromagnetically 
may be described, to terms of second order, by the Darwin Lagrangian |43j . [44j : 



Loarwin = ~<? J2 m i + 3 ^ + fi? 5 ~ 2 S ^ 7^ 

i i i i j^i l 3 

+73 + (vdi^ijVj)) (l) 

where, rrii and are the mass and electric charge of particle i (i = 1, 2, N), Vi its velocity, 
?%i — — fj the relative position of particle i with respect to particle j and c the speed of 

light. In (H}, Hij = ^, 

A generalization of the Darwin Lagrangian that allows the inclusion of other types of 
interactions was proposed by Woodcock and Havas in the 70's [45] . Other formulations have 
recently been studied [46]. For quantum chromodynamics, quark-antiquark interactions to 
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terms of second order have been studied by several authors [17J -[53J- Gravitational interac- 
tions, to terms of second order, are described by the Einstein-Infeld-Hoffmann Lagrangian 



To terms of fourth order (%), ignoring the radiation effects (which for Faraday-Maxwell 
electrodynamics appear at third order) the Lagrangian for an isolated system of N particles 
interacting electromagnetically was obtained in [55], [56]. For General Relativity, gravita- 
tional radiation effects occur at fifth order. To terms of fourth order, the Lagrangian for an 
isolated system of N particles with gravitational interactions was derived in [57]. In recent 
years, for the gravitational two-body problem, results to sixth order (3rd post-Newtonian 
approximation) and higher have been obtained [55] . 

The Lagrangian (CQ) depends only on the velocities and the relative positions of the 
particles, which are functions of a single time variable in a given inertial reference frame 
K. The Darwin Lagrangian can be derived either from Faraday-Maxwell's field theory 
of electrodynamics |44j or from the relativistic action-at-a-distance theory of Wheeler and 
Feynman [2]. 

Wheeler- Feynman's fully relativistic action-at-a-distance theory of electrodynamics is 
described by the Fokker action [1] - [7J: 

S = - Y. m * c J dX iQ - g~ S £ e i e i J J dXidXj^jS (p^) (2) 

In (j2J), Aj is a Poincare invariant parameter that labels the events along the world line 
(Aj) of particle i in flat Minkowski spacetime and, 

C* = % (3) 

iij = ( 4 ) 

pij = - z i? (5) 

are scalars under Poincare transformations. 

The metric tensor: rj^ = diag(+l, —1, —1, —1). In (011]), if = -7^. 

The Dirac delta function in ([2]) accounts for the interactions propagating at the speed of 
light forward and backward in time. 

In action-at-a-distance electrodynamics the interactions carry energy and momentum to 
and from the particles and may simulate a field between them. However, in the absence of 
electrically charged particles outside the system, this fictitious field cannot carry energy or 
momentum into or away from the system. If there are no other electrically charged particles 
in the universe, the total energy, linear momentum and angular momentum of a system of 
point particles are conserved (assuming that only electromagnetic interactions take place). 

In [21], a general action functional depending on the four- vector velocities and relative 
positions of the particles was investigated. Assuming that the Poincare invariant parameters 
that label successive events along the world lines of the particles can be identified with the 
particles' proper times in flat spacetime and that all masses are scalars, the most general 
form of the interaction terms in the action was determined explicitly. 

The theory obtained in [31] is described by the following action functional: 

S = - ^miC J dXiQ - — J2 J2 9i9j J J dXidXjS (p y ) [a^ji + fi&j + iQCf) ( 6 ) 
i i j+i 
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In <Q, a, j3 and 7 are constants, and 



lij = (zi(zj - zt)) (7) 

The Poincare invariant parameters Aj are identified with the proper times (xj) of the 
particles (Aj = Si = cr*) [51] : 



(8) 



In the static limit the theory has three components: a linearly rising potential, a Coulomb- 
like interaction and a dynamical component to the Poincare invariant mass. 

There is strong experimental evidence indicating that, for large separations, the inter- 
actions between quarks can be effectively described by a linearly rising potential [ST], [52] . 
Several relativistic generalizations of a linearly rising potential have been studied [2] - [18] . 
From quantum chromodynamics, it has been shown that the quark-antiquark bound states 
are effectively described by a static potential, which is a sum of a linearly rising potential 
and a Coulomb-like interaction [59] . It would be interesting to study the effect of a variable 
scalar mass for quarks, as given by the formula ([101) below or in some other form. 

From the action ([6]), using the variational principle, one can obtain the relativistic equa- 
tions of motion for a system of N interacting particles [31j : 

(9) 
(10) 



In ([9]) fhi is the dynamical mass of particle % {% = 1, 2, N): 

fhi = rm + 7^ gj J dsjS (pij) 



and is the four- vector force acting on particle i (i — 1, 2, N): 



K? = 9l (Ft u z lu + r 



(11) 



where Ff v is an antisymmetric tensor (F^ = —F^) and V^* is a symmetric tensor (I^g 
r^ a ) in flat spacetime. They are given by the expressions [31] : 



d Sj 



j (Pij 



> Zj ) 2^- ^ ( Z^ Zj 



H [( 

+0 - - - 4)) (1 - - *,))) 

+ Or 



Zj ( z i 



7ji 



(12) 
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It should be emphasized that the action functional ([6]) is the most general expression one 
can obtain from the following four assumptions [3T] : 

(1) the action does not depend on the four- vector accelerations or on higher derivatives, 

(2) the interactions propagate at the speed of light in vacuum, 

(3) the Poincare invariant parameters that label succesive events along the world lines of 
the particles can be identified as their proper times in flat spacetime, 

(4) all masses are scalars under Poincare transformations. 

As it can be seen from ([6]), significant constraints on the form of the action functional 
result from these assumptions. The theory depends only on three undetermined constants. 
It reduces to the Wheeler-Feynman theory of electrodynamics when two of the constants (a 
and 7) are assumed to be equal to zero. 

Based on these results, it is natural to ask how the theory described by the action (jSJ) 
relates to physical phenomena that may be observed experimentally. One possibility is to 
look at the theory as a modified theory of electrodynamics and search for effects arising from 
the new interaction terms (assuming that a or 7 are not exactly equal to zero). Another 
possibility is to consider the theory as an incomplete model for the description of strong 
nuclear interactions (since the theory incorporates a linearly rising potential, but does not 
take into account the internal degrees of freedom of the particles). 

The fully relativistic equations of motion (J9]- [TBI admit exact circular solutions for any 
number of particles. For electrodynamics such solutions were obtained in [6]. 

Our purpose now is to obtain, to terms of second order, the Lagrangian, the Hamiltonian 
and the equations of motion, with all dynamical variables depending on a single time t in 
an inertial reference frame K. To this end, it is convenient to rewrite the action functional 
([6]) in the following form: 



s = -J2 mc 2 Jd U (i-|) 2 -|EE m J J dtidtjd (c\u - t 3 f - (f t - r 3 f) f 13 

(14) 



where, 
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From ffl~4"l [T5l) we can obtain the action for an individual particle i, assuming a fixed 
motion of the other particles, as follows: 



S; 



dt 



-m;C 



1- ^ 



c 9i 2^ 9j 



dtjS^it-tj) 2 - (fi - fj) 2 ) Fi, 



(16) 



where in (fTBj) F i3 - depends on t = and on t 3 (j ^ i). 

The Dirac delta function can be expressed as follows [7]: 
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S^it-tjf-in-rj) 2 ) 



1 

2c 



"3 L j 



6( tj - 1 « 



TDret 



TDadv _i_ V ' J J 



(17) 



In f fT7|) . t^ ,s ' ) (s = — , +) are the two roots of the equation: 

<?(t-t j ) 2 -(7 f i {t)-i f j (t j )) 2 = 



and, 



mf = c{t-t 



R° 



■>:s 

adv 

ij 

f>ret 



Ti — fj 



r>ac 



adv 



Ti — r 



=(+) 



t — tf' ■* is the time it takes for a signal to travel forward in time at the speed of light 
from particle j to particle i in K. 

4 — t is the time it takes for a signal to travel backward in time at the speed of light 
from particle j to particle i in K. 

The action functional ( fT6l) for particle i (assuming a fixed motion of the other particles) 
can therefore be rewritten as: 



(18) 

(19) 
(20) 
(21) 
(22) 
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To terms of second order we can write: 
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We find the approximate action for particle z to be (integration by parts allows to obtain 
a formula without explicit dependence on the accelerations): 



S; 



dtLj 



(27) 



where, 



2 . m i V i . . V- / (^ + 7)' 

-miC + — + - g^- + ft 2. ft" (^ r 'i — 



+^Eft 



— ar,ij + 



(28) 



From (128]) we can find the Lagrangian of the whole system of iV interacting particles. It 
is given by the expression: 



L 



rriivf 



VTiiVi 1 



E-* + E ^ + E ^ + n EEftft ( - 
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+zsEE&& 



ViVj) + ar i3 - + 



(/3 + T)' 



7i, 



A)(nijVj) + — (vf +vfj 



The Lagrangian (129]) reduces to the Darwin Lagrangian ([T]) if a = 7 = and (3=1. 
From the Lagrangian (1291 one can obtain the equations of motion to terms of second 
order: 



(29) 



m. 



70i \ - 9j 



9j 



+ \ a + 



a + 



■I.I 



rl 2c 2 



,2 



rti I \ 2c 2 



- [a + 



2c 2 

3(/? + 7)\ (^) 2 



2c 2 



+ f a + A\ VjirgVi) Vi(rijVi-Vj) + Vjifjjvjf 



+ \a + 



T% J 2c 2 



(30) 



The approximate equations ( 1301 can also be derived directly from the exact equations of 
motion P-H3]). 

From the Lagrangian ([291 we can obtain the total energy (E) , the total linear momentum 

(P) and the total angular momentum (L), of a system of N interacting relativistic particles, 
to terms of second order: 



E 



E^+E^ + E^-SEE^-lar*- 



3m,- 1 



(P + i)' 



i j^i 



Ac 2 



E E 9i9j 



i j^i 



-OLTij + 



09 -7)' 



(ujUj) + ( aTij + — + 7 ^ ) {nijVi){fiijVj) + — (v 2 + 
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(31) 



E \ m i 1 + 



5_ , Tffi^ 0i 
2c 2 / 



E 



\ n ij V j) n ij 



(32) 



2c 2 



7 



i j^i 



-ar ii + 



(0-7)' 



r ij ) r ij 



(fi X fj) 



(33) 



In order to develop the hamiltonian formulation of the theory, we find the conjugate 
momenta: 



Pi 



8L 

dvi 



rrii 



2c 2 rriiC 2 



, 9i \ - 



—arij + 



(0-t)' 



(34) 



The Hamiltonian (to terms of second order) for iV interacting particles is given by the 
formula: 



3^2 2 ^ ^ ' ' " '' 

i j^i 



9i9j 



4c 2 rr! miVriA 

l J jit 1 J 

4c 2 ^ W ™f, 



- ar ii + — ; yPiPj) + ar v + — {nijPi)(nijPj) 



(35) 



For iV = 2 in the center of momentum frame (c.m.f.) (pi = — p~2 = P~) the Hamiltonian 
(|35p takes the simple form: 



H 



c 2 (mi + m 2 ) 



p z 



1 1 

mi m2 



(0 + 7)' 



Z_ ( — + — 

8c 2 \mf m| 



9\92 
2c 2 



m 2 ^ " J 



1 \ 7P 
?«2 I r 



9\92 



2c 2 mim2 



-ar H p 



ar + (/j + 7) ) (np) 2 



(36) 



Direct predictions for the energy spectrum can be made from (1361) . 

Comparing our expression for the hamiltonian (1361) with the general expression presented 
in |3H], we can determine the Lorentz structure of the interactions. 

The general formula for the semirelativistic two-particle hamiltonian presented in [49] 
can be written in the c.m.f. as follows: 



H = c 2 (mi + m 2 ) + 



P 



1 1 
— + — 

mi m 2 



P 



1 1 
-^ + 



nxi 



+ W 



(37) 



where W is the sum of vector and scalar interactions: 

W = W v + W s (38) 
Given a potential V for vector exchange, Wy in the c.m.f. is given by the formula: 

W v = V(r) + = (V(r)p 2 - rV'(r)(npf) (39) 

Given a potential S for scalar exchange, Ws in the c.m.f. can be written as follows: 

Ws = S(r) _ Sir) + 4) g " ^ (S(r)p' + rff(r)(^ (40) 

Comparing ([501) with ([571 -110 )1 we conclude that, for the theory presented in this Letter, 
the potentials corresponding to vector and scalar interactions are given by the following 
expressions: 

V = g ig2 (-ar+?) (41) 

S = *m (42) 
r 

The interaction responsible for the variable masses is scalar in nature. The lineraly rising 
potential is purely a vector potential, and so is the Coulomb-like potential. All this can also 
be seen directly from the fully relativistic action (jSJ). 

For N = 2, the approximate equations of mot ion (1501) admit circular orbits solutions with 
the angular frequency u and the distance between the particles r = r\ + r2 obeying the 
pier rela 



generalized Kepler relation: 



2 0102 

U) = - 

[ir 6 



(ar 2 + (3 + 7 

where, 



9l92 f n _ ... 2 , a ( . . 1 



/j,c 2 r 



Zavr + (3 [v +7^ 



(43) 



, = mim * (44) 
rrii + m 2 

v = 7 Z ^ 45 

(mi + m 2 ) z 

For two-body circular orbits the total energy of the system in the c.m.f., to terms of 
second order, is given by the formula: 



E = c 2 {mi+m2) J ! 



9i92 
2r 



-3ar 2 + (3 + 7 + 



9\92 
4c 2 /ir 



(ar 2 + f3) 2 {\ -v)- 2(ar 2 + pftv - 7 2 (1 + v) 



The total angular momentum in the c.m.f, to terms of second order, is: 



gig 2 {ar 2 + (3 + j)/ir 



1 + 



0102 

4c 2 /Ltr 



Aaur 2 - (3 + 27 



(46) 



(47) 



One can use the procedure presented in this paper to obtain results at higher orders in 
the series expansions. 
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At the quantum level one can describe a system of two spinning particles interacting 
electromagnetically, including second order terms, by the Breit equation [5D]- [22] • An 
extended Breit equation may be devised in agreement with the Lagrangian (1251) to account 
for the dynamical contributions to the masses and the effect of a linearly rising potential. 

It may also be possible to obtain a fully relativistic two-body Dirac equation from the 
action ([6]), by extending the results of Crater and Van Alstine [62] to include the additional 
interactions explored in this Letter. 
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